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ABSTRACT 

The expectations for and the measurements of transverse energy flows, single 
particle px spectra, and the rate of forward jets in deep inelastic ep events 
from the HI and ZEUS experiments at HERA are reported and discussed. It is 
shown that together they offer a good chance to establish deviations from the 
DGLAP paradigm. At the present level of limited statistics the measurements 
are compatible with predictions using BFKL resummation and with the color 
dipole model. Models based on DGLAP evolution describe the pr spectra and 
forward jets less well but are not ruled out yet. 

1. Expectations and Observables 

In deep inelastic scattering (DIS) at low x the simple picture of the quark parton 
model, where a virtual photon interacts with a point-like parton in the proton and 
nothing else happens, has to be modified. The probability that additional partons, 
particularly gluons, are radiated increases with decreasing x. An example, where a low 
X parton, which interacts with the photon, originates from a parton shower initiated 
by a gluon with large x, is shown in Fig. |l|. In the approximation of just one gluon 
being radiated in the initial or the final state, such processes have been calculated 
in next to leading order (NLO). Two programs, MEPJET and DISENT i, using 
different methods, are available, and a new program, DISASTER++ i, has been 
announced at this meeting. These programs generate partons only; no full event 
generators exist. 

In the leading log approxima- 
tion any number of more or less 
collinear or soft gluons may be ra- ^ 1 

diated. By resumming terms of the si 
form (as In C)" using the DGLAP § 

evolution equations an initial state s! '^i' '^t \ 

parton shower may be evolved from 
the scale Ql, typically between 0.3 
and 4.0 GeV^, to the scale Q^, 
where the interaction with the pho- 
ton takes place. DGLAP evolution 
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Fig. 1. Diagram of an ep collision at low x. 
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implies weakly ordered fractional gluon energies (longitudinal momenta) Xi > Xi+i 
and strongly ordered gluon virtualities or transverse gluon momenta fcy,* << ^t.j+i 
along the ladder (Fig. [1|). Different implementations of "DGLAP-like" parton showers 
and of matching EeI them to leading order (LO) exact matrix elements for the QCD- 
Compton (QCDC) and the boson-gluon fusion (BGF) process exist in LEPTO 0, 
RAPGAP i, and HERWIG i. For the hadronisation of the perturbative partons to 
the observable hadrons, these programs include phenomenological models. HERWIG 
uses the cluster model 0; all other programs, including ARIADNE 0, use the LUND 
string model lil, as implemented in the JETSET El code. These fragmentation models 
have been well tested at LEP 0. At HERA and at hadron colliders, the additional 
complication of proton remnant fragmentation arises, about which much less is known. 

At low X, terms of the form (a^ln^)" become large and can be resummed using 
the BFKL Il3 evolution equations. In a physical gauge, these terms correspond to a 
ladder diagram in which gluon emissions are strongly ordered in energies Xi >> Xi+i, 
while, in contrast to the DGLAP case, the transverse gluon momenta perform a 
random walk along the ladder. Analytical results at the parton level and first results 
including fragmentation functions exist and will be mentioned later. Unfortunately 
no implementations of "BFKL-like" parton showers are available in any of the event 
generators. 

Finally we discuss ARIADNE 0, which currently gives the best overall description 
of the hadronic final state in DIS, but which does not fit well into anv of the schemes 
mentioned above. It is based on the color dipole model (CDM) EJ, where gluon 
emission originates from a color dipole stretched between the scattered quark and the 
proton remnant. Further gluon emission leads to a chain of independently radiating 
dipoles spanned between color-connected partons. The first emission is corrected by 
the QCDC matrix element in LO. In ARIADNE, the color charge of the proton 
remnant is assumed not to be point-like, leading to a phenomenological suppression 
of gluon radiation in the remnant direction. The suppression sets in for hard 
gluons with wavelengths smaller than the size of the remnant. Also the color charge 
of the scattered quark is taken to be spread out, depending on the virtuality of 
the photon 00. This in turn leads to a suppression of radiation in the direction of 
the scattered quark. An important feature of ARIADNE is that the probability for a 
QCDC event and for any other gluon emission, from a dipole connected to the remnant 
at one end, depends on this "unorthodox suppression" 0, due to the extended color 
charge, while in "normal" QCD, as implemented in all other programs, it depends 
on the ratio of parton densities before and after the emission and on other factors 
which are the same for all programs. At low x this feature leads to a greatly enhanced 
rate of QCDC events compared to the other generators. Additionally, in contrast to 
the DGLAP-like programs, the gluons in ARIADNE are not ordered in kf along the 
ladder. It has therefore been argued that the partonic final state of ARIADNE is 
more closely related to the one expected from BFKL evolution. 



It has been shown that all models provide a fair description of basic event prop- 
erties Bi'i. 

At very low x the gluon density may possibly become so high in regions ("hot 
spots") of the proton that the gluons will no longer act as free particles but will 
interact with each other. In this saturation regime of QCD, the strong coupling is 
still small, but the incoherent scattering approximation is no longer valid. It is not 
clear whether this new regime could be observed at HERA. 

What do we expect to see at HERA? Which observables will be most sensitive to 
the underlying parton dynamics? 

As X decreases below 10~^ we anticipate deviations from the DGLAP predictions 
due to the missing (agln^)" terms. The most inclusive way to study the structure 
of the proton is to measure the probing lepton after the interaction. The resulting 
measurement of the proton structure function F2 has the advantages that it is ex- 
perimentally relatively easy and that it can be directly compared to analytical QCD 
calculations. The rise of F2 with decreasing x as observed at HERA §1 turns out to be 
well describable by the DGLAP evolution equations alone, down to > 1.0 GeV^ 
and X > 10~^. At this meeting, it was shown that the F2 data at small x are also 
well compatible with a unified BFKL and DGLAP description 0. 

Perhaps F2 is not sufficiently sensitive to observe the change in parton dynamics 
at low X. The reasons might be that F2 is a too inclusive measurement, averaging 
out over differences only visible in the details of the hadronic final state. There might 
not be enough range in x at moderate Q^. Moreover, the theory has some fiexibility 
due to the freedom of choice of the starting parton distributions. 

A greater sensitivity to small x effects might be provided by observables based 
on the hadronic final state emerging from the underlying partons. It can be further 
enhanced by making use of the differences between the DGLAP and the BFKL cas- 
cade. However, there is a price to pay. The measurements have to be compared to 
Monte Carlo simulations including the QCD effects, the transition from partons to 
hadrons, the particle decays, and detector effects. The observables studied so far are 
the transverse energy flows, single particle pT-spectra, and forward jets. 

2. Transverse Energy Flows 

BFKL dynamics predicts HIl, due to the lack of ordering of kr in the parton 
cascade, more average transverse energy Ft and harder tails in the Ft distribution 
for decreasing x than in the DGLAP scenario. These observables are measured calori- 
metrically in the central rapidity " region, between the struck quark and the proton 
remnant, in the hadronic center of mass (hems), that is the rest system of the ex- 
changed virtual photon and proton. In that frame all Ft * is due to either QCD 

"Rapidity is here always pseudo-rapidity, i.e. rj = — lntan(0/2). 
''Et = ^Ei sin. 6i, i runs over calorimeter cells. 



radiation or hadronisation. As indicated in Fig. |1| the central rapidity in the hems 
corresponds to the forward region in the frame of the HERA detectors. 




Fig. 2. Transverse energy flow as a function of rapidity 77* in the hems. The HI data points show 
statistical errors only, except for two points measured with the plug calorimeter where the statistical 
and systematic errors are given. The lines are the Monte Carlo predictions of three different models. 

First results on the Ex-fiows have been published by the HI collaboration BSil. 
In ref. ^ they have been reported for DIS events with 5 GeV^ < < 50 GeV^ 
and with 10~^ < x < 10~^ in the hems and have been compared to the then current 
model predictions of LEPTO 6.1 and ARIADNE 4.03. In the largest x and bin, 
i.e. x ~ 5 X 10~^ and ~ 33 GeV^, the data and the two models were found to 
agree. With decreasing x for fixed the Ex from LEPTO 6.1 was found to fall more 
and more below the data for rapidities away from the current jet, while ARIADNE 
4.03 still managed to give a level of E^ in agreement with the data. 

At this point in time one could have had the impression that models with DGLAP- 
like dynamics like LEPTO fail to describe the data at low x, while models with some 



BFKL-like features like ARIADNE are successfull. However, soon after, two new 
phenomenological features have been added to LEPTO which allowed a reasonably 
good description of the transverse energy flows also for decreasing x. The new features 
are soft color interactions S, with the intention to describe rapidity gap events without 
modeling a Pomeron and its structure function, and a modified sea-quark/remnant 
treatment, giving a smoother transition from BGF events to events where the photon 
interacts with quark. 

The HI collaboration in the meantime has produced new preliminary data 011, 
covering a larger range in the kinematic plane from 3 < < 70 GeV^ and 8 x 10~^ < 
X < 7 X 10~^. The extension to lower was achieved by analysing data from special 
runs, where the point of the ep interactions was shifted from the nominal position in 
the direction of the proton beam in order to have access to smaller lepton scattering 
angles. In addition, the Et at two very forward rapidity bins was measured by HI 
using their plug calorimeter (0.72° < 9iab < 3.3°), which closes the gap between the 
beam-pipe and the forward part of their liquid Argon calorimeter. 

The data as a function of rapidity in the hems are shown in bins of x and in 
Fig. 0. They are compared to predictions from ARIADNE 4.08, HERWIG 5.8, and 
LEPTO 6.4. The models give only a fair description of the data over the large range in 
X and Q^. The average Et in the central rapiditv —0.5 <•)]*< 0.5 as a function of x in 
bins of fixed increases with decreasing x eIeHI. This is demonstrated in Fig. ^ for 
the bin = 14 GeV^, together with data from ZEUS il, and with predictions from 
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Fig. 3. Mean transverse energy in the central rapidity bin in the hems for the bin = 14 GeV^. 
Besides the data from HI and ZEUS, QCD model predictions (lines) are shown for hadrons in a) 
and for partons in b). An analytic BFKL calculation at the parton level is shown as open circles. 



different generators at the hadron level. For LEPTO the predictions are shown with 
the new soft color interactions (SCI) and the new sea-quark treatment (SQT) turned 
on and off. The data of both experiments are in good agreement and can be described 



■^The direction of the proton is to the left (negative rapidity). 



by all models, DGLAP and BFKL-like. In Fig. ^ the predictions of the models at the 
hadron level are contrasted with those at the parton level and with an analytic parton 
level BFKL calculation @. At the parton level the DGLAP-like models LEPTO and 
HERWIG show the opposite slope in x than ARIADNE, the BFKL calculation, the 
hadron level of all models, and the data. The BFKL calculation gives the highest 
transverse energy, 40% to 50% above the partons from the color dipole model. There 
are large differences in the contribution from fragmentation to the total Et between 
the models. It is about 15% for the BFKL calculation (comparing it to data), 35% 
for ARIADNE, and about 70% for LEPTO and HERWIG. 




E-r (GeV) 



Fig. 4. Transverse energy distribution in the central rapidity bin in the hems. Only statistical errors 
are shown. 

Concerning the measurement of the mean Et in the central rapidity bin, we are 
now in the situation that all models describe the data. Although the models differ in 
their underlying parton dynamics, the mean E^ can be made to agree by exploiting 
as yet unconstrained variations in hadronisation models. 

Another observable investigated by HI is the distribution of the Et in the central 
rapidity bin per event. The Et originating from hadronisation is expected to be 
limited, while high values of Et are more likely to be produced by hard parton 



radiation. Preliminary data llB, in the same bins of x and as in Fig. ^, are 
displayed in Fig. | together with predictions from ARIADNE 4.08, HERWIG 5.8, 
and LEPTO 6.4. While at the largest x and the data and the models agree, with 
decreasing x and the data and ARIADNE appear to exhibit harder tails than the 
DGLAP-like models. As can be seen from the figure a firm conclusion can only be 
drawn with more statistics extending the range in Ex to larger values. Comparing 
the Et distribution for the highest and lowest x value for fixed Q^, one finds that the 
data have harder tails with decreasing x, suggesting a rise in parton activity. 



3. Single Particle px Spectra 



The Et measured by the calorimeter can be due to many soft particles, predom- 
inantly from hadronisation, or due to a few hard particles, mainly from hard gluon 
emission, or due to both. Therefore it has been suggested that the hard tail of the 
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Fig. 5. The transverse momentum spectra of charged particles in the hems. The data are shown for 
nine different kinematic bins and the combined sample (bin 0). Statistical and full errors are given. 



Pt distribution of single charged particles in the hems might offer better sensitivity 
to the basic parton dynamics since it is more directly linked to hard gluon emissions. 

The charged particle pt distributions in bins of x and as measured by HI 
in the rapidity interval from 0.5 to 1.5 are shown in Fig. ^. At large x all three 
models presented agree with the data. With decreasing LEPTO and HERWIG 
significantly fall below the data for increasing px- Predictions based on BFKL 
resummation and convolution with fragmentation functions for the transition from 
partons to charged hadrons § have been made for the three lowest bins in i.e. 
bins 1, 2, and 3. As demonstrated in Fig. ^ for bin 3, the data and the BFKL 
result agree quite well for px > 1.5 GeV. The absolute normalization was derived 
from a comparison of the BFKL calculation of the forward jet cross section with data 
from HI to be discussed later. Also displayed in the figure is the calculation with 
BFKL effects turned off, which falls below the data. It is apparent from the figure 
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Fig. 6. The transverse momentum spectra of charged particles in the hems for fixed low x and 
(bin 3). The data and predictions from a BFKL calculation, with BFKL effects off and on, and with 
two different choices for the scale /i of the fragmentation function are compared. 



that the significance of the agreement can be increased with more data at higher 
Pt- In addition, BFKL effects would become stronger, if the measurement of charged 
particles could be extended further in the direction of the proton remnant, i.e. to the 
rapidity interval —0.5 < r]* < 0.5. 



4. Forward Jets 



The cross section for forward jets in DIS as a function of x has been advertised £3 
for some time now as an observable enhancing the effects of BFKL resummation. 
Diagrammatically the situation is described in Fig. 0. The forward jet is defined by 



the azimuthal angle 6jet between the jet and the proton direction and its transverse 
momentum fcrjet- The energy of the jet is then given 
by Ejet = kTjet/ sin 6 jet and xjet ~ Ejet/Ep. We now re- 
quire that k^jet ~ which suppresses the phase space 
for forward jet production and for gluon emission be- 
tween the forward and the current jet in the DGLAP 
case. In addition, we demand that xjet/x he as large as 
possible which maximizes the phase space for the pro- 
duction of a forward jet and gluon radiation between 
it and the current jet in the BFKL scenario. For this 
particular kinematic configuration the resummation of 
as In Xjet/x terms should lead to a sizeable growth of the 
forward jet cross section with decreasing x. Another ad- 
vantage of this observable is that the parton densities 
of the proton are probed at rather large scales, xjet and 
Q^, where they are well known. 
Published results on forward jets are available from HI H and new preliminary 
results from HI H and from ZEUS @. The cuts for the selection of DIS events and the 
requirements on the forward jet are summarized in Table |l]. Both experiments use 
a cone algorithm, requiring a minimum prjet in a cone of radius R = \/ Ai]"^ + A$2 
in the HERA frame. There are differences in the current analyses between the two 




Fig. 7. Parton evolution with 
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Table 1. DIS event and forward jet selection 





HI 


ZEUS 


Ee > 


11 GeV 


10 GeV 


Ve > 


0.1 


0.1 


0e< 


173° 




Oe > 


160° 




cone in lab, R 


1.0 


1.0 


djet > 


7° 


10° 


PTjet > 


3.5 GeV 


5.0 GeV 


Xjet > 


0.035 


0.035 


Sjet < 


20° 




PTjet ^ 


0.5 


0.5 


PTjet ^ 


2.0 


4.0 



experiments in the cuts on Qj^t-, PTjet-, and the upper limit on p^jetlQ^ which will 
be discussed later. With the CDM the transverse energy flow around the forward 
jet axis can be well described for different values of PTjet ^ as well as many other 
distributions 



In Fig. 1^ the ZEUS forward jet cross section corrected to the parton level of 
ARIADNE is shown and compared to several parton level calculations: an analytic 
BFKL calculation 0, the same calculation but without any gluon emission between 
the forward and the current jet system (Born BFKL), and a fixed NLO QCD cal- 
culation using MEPJET. The data show a much faster rise with decreasing x than 
the calculations without BFKL resummation. The BFKL calculation shows an even 
more dramatic rise. The authors 0, however, point out that several effects which 
have not been taken into account might lower their prediction. ZEUS presents their 
preliminary result corrected to the parton level in order to compare to parton level 
calculations. Using ARIADNE correction factors are found to vary between 0.6 for 
the lowest x bin and ~ 1 for larger x. It will be useful to have the ZEUS cross section 
also at the hadron level, since it is not clear, what the relationship is between the 
parton level of the CDM and BFKL. It could also allow comparisons with HI data 
which are corrected to the hadron level. 
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Fig. 8. The difFerrential forward jet cross section as a function of x. The ZEUS data have been 
corrected to the parton level. Statistical and systematic errors (not yet complete) are included. 
BFKL calculations, with BFKL effects on (dash-dotted line) and off (dotted line), calculations in 
NLO QCD using MEPJET, and parton level results from the color dipole model are shown. 



The preliminary HI results on the forward cross section are shown in Fig. ^ 
and are compared to predictions from LEPTO (MEPS) with and without soft color 
interactions (SCI) and to ARIADNE (CDM). Again a fast rise of the cross section 
is observed with the CDM falling only slightly below the data. In Fig. models 
and calculations at the parton level are shown: ARIADNE, LEPTO, two BFKL 
calculation and a NLO calculation using DISENT. ARIADNE shows a similar 
X dependence at the parton and hadron level with hadronisation effects amounting 
to less than 20%. The parton level forward jet cross sections of LEPTO and of 



DISENT in NLO agree as expected and show a moderate increase. In LEPTO, 
with soft color interactions turned on, up to 80% of the forward jets are created in 
the hadronisation phase causing LEPTO to only slightly undershoot the data with 
decreasing x. The BFKL calculations basically can only predict the dependence on 
X but not the normalization. As mentioned before, the computation by Kwiecinski 
et al. @ fixed the normalization to the data from HI while the calculation by Bartels 
et al. i did not. Scaling the cross section by Bartels et al. down by a factor 0.8 
brings the two BFKL calculations in agreement. The color dipole model provides the 
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Fig. 9. The forward jet cross section as a function of x. The HI data have been corrected to the 
hadron level. Statistical and systematic errors are included. QCD model predictions for hadrons 
are superimposed in (a) and for partons in (b). Two analytic BFKL calculation (open squares and 
triangles) and a NLO calculation using DISENT (open circles) are also shown in (b). 

best description of the data. However, it should be pointed out that its prediction is 
rather sensitive to the power of soft suppression of gluon emission {^/pt)" due to the 
proton remnant. The parameter a is related to the dimensionality of the extended 
proton remnant and is expected to be ~ 1. The HI forward jets prefer a < 1.0 while 
non-jet observables like energy flows prefer a ~ 1.5 0. 

One may wonder about the choice for the minimum pTjet, 3.5 GeV by HI and 5.0 
GeV by ZEUS. A small value is desirable from the point of view of statistics and 
sensitivity to BFKL effects @. With increasing prjet the slope in x for the forward jet 
cross section decreases. The requirement prjet ~ forces the to increase, which 
in turn increases the minimum x which is probed. On the other hand, to suppress 
hadronisation effects, one would want prjet to be large. 

Another question concerns the choice of the upper limit on Pj'jet/Q'^ which is 2.0 
for HI and 4.0 for ZEUS. Increasing this limit from 2.0 to 4.0 increases in the lowest 



X bins the contribution from 0{as) matrix elements to the forward jets by roughly 
a factor of two. One would also expect an increased contribution from resolved 
photoproduction which is not included in the calculations. After this workshop the 
author of RAPGAP has shown S that for example the HI result can also be described 
with direct (same as in LEPTO without SCI) and resolved contributions. 

5. Conclusions 

The tails of the single particle pt spectra and of the calorimetric Et distribution, 
and the forward jet cross sections offer a good chance to pin down deviations from the 
DGLAP paradigm. Other interesting observables in this quest, like the decorrelation 
of the azimuthal angle between the forward jet and the leptonc^, the forward tc° cross 
section 0, and the production of more than one forward jet should be pursued. 

Along the way, we probably will and have to get a better understanding of hadro- 
nisation effects, particularly of the proton remnant. On the theory side, higher order 
effects @ have to be included in the calculations and a BFKL Monte Carlo genera- 
tor iiS is needed. 

In addition to collecting more data at HERA, it would be desirable to access 
smaller forward angles for jets E and particles and to increase the HERA center of 
mass energy. 

6. Acknowledgements 

I want to thank T. Carli, M. Wobisch, and S. Wolfle for providing me with plots 
and J. Dainton, J. Hartmann, and D. Kriicker for a careful reading of the manuscript. 
I am also happy to thank the organizing committee, in particular B. Kniel, for the 
excellent workshop and the pleasant atmosphere. 

7. References 



1. E. Mirkes and D. Zeppenfeld, Phys. Lett. B380 (1996) 205. 

2. S. Catani and M. Seymour, Proc. Workshop on Future Physics at HERA 
1995/96, Hamburg, eds. G. Ingelman, A. De Roeck, and R. Klanner, vol. 
1 (1996) 519, CERN-TH/96-240 (1996) [|hep-ph/9609521|| , and Acta Phys. 
Polon. B28 (1997) 863. 

3. D. Graudenz, these proceedings. 

4. V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438 and 675; 
G. Altarelh and G. Parisi, Nucl. Phys. B126 (1977) 298; 

Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

5. M. Bengtsson, G. Ingelman, and T. Sjostrand, Proc. HERA Workshop 1987, 
Hamburg, ed. R.D. Peccei, vol. 1 (1988) 149; 



M. Bengtsson and T. Sjostrand, Z. Phys. C37 (1988) 465. 

6. M. H. Seymour, LU-TP-94-12, Nucl. Phys. B436 (1995) 443 and Comput. 
Phys. Commun. 90 (1995) 95. 

7. G. Ingelman, Proc. Workshop on Physics at HERA 1991, Hamburg, eds. W. 
Buchmiiller and G. Ingelman, vol. 3 (1992) 1366; 

G. Ingelman, A. Edin, and J. Rathsman, Comput. Phys. Commun. 101 (1997) 
108. 

8. H. Jung, Comput. Phys. Commun. 86 (1995) 147. 

9. G. Marchesini et al., Comput. Phys. Commun. 67 (1992) 465. 

10. B. R. Webber, Nucl. Phys. B238 (1984) 492; 

G. Marchesini and B. R. Webber, Nucl. Phys. B310 (1988) 461. 

11. L. L5nnblad, Comput. Phys. Comm. 71 (1992) 15 and Z. Phys. C65 (1995) 
285. 

12. B. Andersson, G. Gustafson, G. Ingelman, and T. Sjostrand, Phys. Rep. 97 
(1983) 31. 

13. T. Sjostrand, Comput. Phys. Commun. 82 (1994) 74; 

T. Sjostrand, PYTHIA 5.7 and JETSET 7.4, CERN-TH 7112/93 (1993) and 
LU TP 95-20 (1995). 

14. I.e. Knowles et al., Proc. LEP2 Physics Workshop 1995 ||hep-ph/ 9601211 . 

15. E. A. Kuraev, L.N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 45 (1977) 199; 
Ya. Ya. Balitsky and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822. 

16. G. Gustafson, Phys. Lett. B175 (1986) 453; 

G. Gustafson and U. Petterson, Nucl. Phys. B306 (1988); 

B. Andersson, G. Gustafson, L. Lonnblad, and U. Petterson, Z. Phys. C43 

(1989) 625. 

17. L. Lonnblad, M. Seymour, et al., Proc. Physics at LEP 2 Workshop, eds. G. 
Altarelh, T. Sjostrand, and F. Zwirner, CERN/96-01, vol. 2, (1996) 187. 

18. J. Rathsman, Phys. Lett. B393 (1997) 181. 

19. A.H. Mueller, Nucl. Phys. B415 (1994) 373. 

20. N. Brook et al., Proc. Workshop on Future Physics at HERA 1995/96, Ham- 
burg, eds. G. Ingelman, A. De Roeck, and R. Klanner, vol. 1 (1996) 613. 

21. T. Carli, Proc. 4th Int. Workshop on Deep Inelastic Scattering and Related 
Phenomena (DIS 96), Roma, eds. G. D'Agostini and A. Nigro, (1997) 415. 

22. H. Jung, private communication. 

23. L. Favart, these proceedings. 

24. A. D. Martin, these proceedings. 

25. K. Golec-Biernat, J. Kwiecinski, A. D. Martin, and P. J. Sutton, Phys. Rev. 
D50 (1994) 217 and Phys. Lett. B335 (1994) 220. 

26. J. Bartels, H. Lotter, and M. Vogt, Phys. Lett. B373 (1996) 215. 

27. HI Collaboration, T. Ahmed et al., Phys. Lett. B298 (1993) 469. 

28. HI Collaboration, I. Abt et al., Z. Phys. C63 (1994) 377. 



29. HI Collaboration, S. Aid et al, Pliys. Lett. B356 (1995) 118. 

30. A. Edin, G. Ingelman, and J. Rathsman, Z. Phys. C75 (1997) 57 and Phys. 
Lett. B366 (1996) 371. 

31. HI Collaboration, C. Adloff et al, contrib. paper pa02-073, ICHEP '96, War- 
saw, Poland, July 1996. 

32. M. F. Hess, Ph.D. Thesis, University of Hamburg (1996). 

33. N. A. Pavel, Proc. 4th Int. Workshop on Deep Inelastic Scattering and Related 
Phenomena (DIS96), Roma, eds. G. D'Agostini and A. Nigro, (1997) 502. 

34. M. Kuhlen, Phys. Lett. B382 (1996) 441. 

35. HI Collaboration, C. Adloff et al., Nucl. Phys. B485 (1997) 3. 

36. J. Kwiecinski, S. C. Lang, and A.D. Martin, DTP-97-56 (1997) [ |h^ 
ph/9707240| . 



37. J. Binnewies, B. A. Kniehl, and G. Kramer, Phys. Rev. D52 (1995) 4947. 

38. A.H. Mueller, Nucl. Phys. B (Proc. Suppl.) 18C (1990) 125 and J. Phys. G17 
(1991) 1443; 

W. K. Tang, Phys. Lett. B278 (1992) 363; 

J. Bartels, A. De Roeck, and M. Loewe, Z. Phys. C54 (1992) 635; 

A. De Roeck, Nucl. Phys. B (Proc. Suppl.) 29A (1992) 61; 

J. Kwiecinski, A. D. Martin, and P. J. Sutton, Phys. Lett. B287 (1992) 254, 

Nucl. Phys. B (Proc. Suppl.) 29A (1992) 67, and Phys. Rev. D46 (1992) 921. 

39. J. G. Contreras, Proc. XXXI Rencontres de Moriond in QCD and High Energy 
Hadronic Interactions, ed. J. Tran Thanh Van, Edition Frontieres (1996) 411; 
HI Collaboration, C. Adloff et al., contrib. paper pa03-049, ICHEP '96, War- 
saw, Poland, July 1996. 

40. S. Wolfie, to be published in Proc. 5th Int. Workshop on Deep Inelastic Scat- 
tering and QCD, Chicago, April 1997, eds. D. Krakauer and J. Respond. 

41. J. G. Contreras, Ph.D. Thesis, University of Dortmund (1997). 

42. E. M. Lobodzinska, Ph.D. Thesis, Inst, of Nucl. Physics, Cracow (1997). 

43. J. Bartels et al, Phys. Lett. B384 (1996) 300. 

44. J. Bartels, A. De Roeck, and M. Wiisthoff, Proc. Workshop on Future Physics 
at HERA 1995/96, Hamburg, eds. G. Ingelman, A. De Roeck, and R. Klanner, 
vol. 1 (1996) 598. 

45. H. Jung, to be published in Proc. Madrid Workshop on Low x Physics, Mi- 
raflores de la Sierra, June 1997. 

46. J. Kwiecinski, S. C. Lang, and A.D. Martin, Phys. Rev. D55 (1997) 1273. 

47. J. Kwiecinski, C. A. M. Lewis, and A.D. Martin, DTP-97-58 (1997) [|Eep 
ph/9707375| 



48. J. Bartels, these proceedings. 

49. C. R. Schmidt, Phys. Rev. Lett. 78 (1997) 4531. 

50. L.H. Orr and W. J. Stirling, DTP-97-48 (1997) ||hep-ph/9706529| , to be pub- 
lished in Phys. Rev. D. 



